Correlation Inequalities for Quantum Spin Systems with Quenched Centered Disorder 
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It is shown that random quantum spin systems with centered disorder satisfy correlation inequali- 
ties previously proved [l]] in the classical case. Consequences include monotone approach of pressure 
and ground state energy to the thermodynamic limit. Signs and bounds on the surface pressures 
for different boundary conditions are also derived for finite range potentials. 
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Quantum spin system with quenched randomness are 
both important and theoretically challenging. They are 
widely used as models for metallic alloys in condensed 
matter physics, see for a review. They are also impor- 
tant in combinatorial optimization problems especially in 
relation to quantum annealing procedures [3] and quan- 
tum error correcting codes f3| . 

One approach to studying such macroscopic systems 
is via correlation inequalities. These are useful in many 
areas of statistical mechanics. They have been used to 
prove that the free energy, correlation functions and sur- 
face free energy have well defined thermodynamic lim- 
its. They further show that some of these quantities ap- 
proach their limit monotonically at increasing volumes. 
They have also proved important for computing bounds 
on critical temperatures and critical exponents by com- 
paring lattices in different dimensions. 

In a previous work [1] we have shown that, despite 
the presence of competing interactions, a general classi- 
cal spin glass model with quenched centered disorder has 
a family of positive correlation functions. This positivity 
implies a monotone behavior of the pressure with respect 
to the strength of the random interaction. From the 
monotonicity one can deduce sub-additivity of the free 
energy (which implies existence of its thermodynamic 
limit) and also the sign of and bounds on the surface 
pressure for different boundary conditions. See also Q 
for further applications of the inequalities. 

The extension to the quantum case of classical corre- 
lation inequalities may not be possible or require further 
conditions. Examples were found in Q where some cor- 
relations violate the GKS inequality of type II for the 
isotropic Heisenberg model with ferromagnetic interac- 
tions, and in [lo| where the GKS inequality of type I 
is violated for the anisotropic ferromagnetic Heisenberg 
model. In order to reestablish the validity of the GKS 
inequalities in quantum systems it is then necessary to 
impose further conditions on the interaction coefficients 
beyond positivity (see for instance [lol | for conditions to 
prove GKS I in the anisoptropic Heisenberg model with 
zero magnetic field and for the conditions to prove 



GKS II). 

Here we show that quantum systems with quenched 
centered disorder do fulfill the same family of correlation 
inequalities of the classical case without any further re- 
striction with respect to the classical case. The result is 
obtained as follows. 

For each finite set of points A let us consider the quan- 
tum spin system with Hamiltonian 



U 



XxJx'i'x + Uo 



(1) 
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The operators are self-adjoint elements of the real 
algebra generated by the set of spin operators, the Pauli 
matrices, a^^\a^^\al^\ i G A, on the Hilbert space 
Tix '■= (Sii^x'Hi. Uq is a non random quantum Hamil- 
tonian acting on the Hilbert space Ha. The random in- 
teractions Jx are centered and mutually independent i.e. 
Av{Jx) = for all X and Av{JxJy) = '!^xSx,y- The 
A's are numbers which tune the magnitude of the random 
interactions. An example is the anisotropic quantum ver- 
sion of the nearest-neighboor Edwards-Anderson model 
with transverse field. This is defined in terms of the Pauli 
matrices: 



(2) 



(3) 



for |« — j| = 1 and ^x = otherwise. 

Our main observation is that the pressure (Gibbs free 
energy up to a sign) for A = (Ax, Ay, ...): 



Pa(A) = AwlogTr exp(-C/), 



(4) 



is convex with respect to each Xx- We have set the in- 
verse temperature /3 = 1 since our results do not depend 
on its value. We shall also drop the subscript A when it 
is unambiguous. 

The proof of convexity is straightforward. The first 
derivative gives in fact 

dP 

— — = Av{Ja < $A >;7) (5) 
oXa 
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where 



<C>[/:= 



TrCe- 
Tr e~' 



(6) 



while, for the second derivative, one has (see Chap- 
ter IV, page 357) 



Av{Jl[< $A, <^>A>u - < *A >?/]) (7) 



where < •, • >[/ denotes the Duhamel inner product [13| : 

<C,i?>a:=— ^r^Ti^ • (8) 

By using the fact that < 1,D >=< D > and < C, 1 >= 
< C > we see that 

02 p 

Av{jl[< $A- < <^A >U,'^A~ < >U>U]) > 0. 

(9) 

This yields the following result: 

For systems described by the quantum potential ([T]) the 
following inequality holds: for all A d A and for Xa > 



Av{Ja <<i>A>u)>0. 



(10) 



Proof. Since the second derivative of the pressure is non 
negative 



> 0. 



we deduce that the first derivative 
dP 

Av{Ja < <^A >u) 



OXa 



(11) 



(12) 



is a monotone non decreasing function of Xa (indipcn- 
dently of the values of all the other A's). As a conse- 
quence we have that for > 



dP 
dX^ 



> Av{Ja < $a > 



U)\Xa=0 



(13) 



But for Xa = the two random variables Ja and < 
^A >u are independent: 

Av{Ja <^a> \xA=a) = Av{Ja)Av{< $a > |a^=o) = 

(14) 

where the last equality comes from having chosen distri- 
butions with Av{Ja) — 0. It also follows that for < 
one has Av{Ja < $a >u) < 0. 

Although the consequences we are going to derive ap- 
ply only to the case considered in ([T]) where Uq is the 
sum of one body terms we note that the inequality 
holds for general Uq. This include the case where 
Av{< ^A >u) £ 0, as would happen in the case where 



the Jx are bounded and U satisfy the conditions neces- 
sary for GKS I to hold. A different example where one 
exploits symmetry and translation invariance would be 
the anisotropic Heisenberg model 



U = - 



(h 



(15) 



ct—x^y^z 



with centered Ji and negative field h. It would also in- 
clude the case /i = 0, A/'Z'^,d>3 with minus bound- 
ary conditions and Ka positive and large. 

We now consider the case where Uq is a sum of one 
body terms, e.g. Uq = — hi ■ Oi. By using the same 
standard strategies of the classical spi n gl ass case [H or 
the standard ferromagnetic interaction [12l | one can easily 
deduce from pO)) the super-additivity of the pressure. 
For a disjoint union of two regions A = Ai U A2 one 
obtains 



Pa > Pai + Pa. 



(16) 



It follows from (|16p that the pressure is monotonically in- 
creasing as the volume increase and hence the existence 
of the thermodynamic limit (see also [3])- Considering 
for instance a system on a d-dimensional square lattice 
Z'', with translation invariant distributions of the ran- 
dom interactions, one has that by dividing the lattice 
into cubes the following result holds for free boundary 
conditions: 



V 



= Hm -j-^ 

A/'Z'i |A| 



Pa 



(17) 



where the supremum is a well defined function (doesn't 
blow up) provided the stability condition (see 14|): 



E 
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Av{Jl)\\^xf < c|A|, 



(18) 



is verified for some positive constant c. A simple bound 
shows that when the interactions have a finite range the 
limit does not depend on boundary conditions. 

By introducing the inverse temperature in the defini- 
tion of the pressure, for instance taking all the lambdas 
equal to /3, we can study the properties of the ground 
state energy E/^ by relating it to the free energy. Since 
by general thermodynamic arguments (see for instance 



one obtains 



which implies 



lim — 

0^00 



Ea 



Ea < Ea, + Ea, , 



Ea . ^Ea 
e = lim — = mt — 

A/'Z'' |A| A |A| 



(19) 



(20) 



(21) 
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The physical significance of a quantum disordered model 
is related to the fact that the random free energy (using 
now the potential with all \x ~ 1) 



nA-logTrexp(-C/A) 
and the random ground state energy 



Ejy = hm — — 



(22) 



(23) 



do converge, for large volumes, to the same non random 
object for almost all the disorder realizations. Following 
[IJ and [17| we can achieve this stronger version of the 
existence of the thermodynamic limit by observing that 
the condition (fT8|) entails the exponential version of the 
law of large numbers for the free and ground state energy: 



Prob 



Prob 



Ha Pa 



IMP |A|/3 

£a _ -Ba 
|A| |A| 



> x] <e~ 



> x] < 



(24) 



(25) 



Standard probability theory (Borel-CantcUi lemma) im- 
plies that for almost all configurations of the J's: 



hm — 

A/'Z'^ \A\ 



Pa 



P : 



and 



Sa . ,Ea 

hm TTT = mi TTT = 

Ayz-i |A| A |A| 



(26) 



(27) 



From inequality (jlOp we can also deduce, by decompos- 
ing a (d-l-l)-dimensional hypercube A into d-dimcnsional 
hypercubes, that the pressure in dimension d, p^"^^ is a 
non-decreasing function of d and the ground state energy 
a non increasing one: 



P 



(d) < p{d+l) 



(28) 



(29) 



In the case of finite range interactions, i.e. Ax = for 
\X\ > r (e.g. the nearest neighbor case), the inequality 
leads to an estimate of the size and sign of the surface 
pressures Ta i.e. the first correction to the leading term 
of the pressure: 



Pa =p|A|-l-rA 



(30) 



Using the methods of [H and [l^ a straightforward com- 
putation shows that the Ta is of surface size and, as it 
happen for ferromagnets and classical spin glasses, it does 
depend on boundary conditions. For instance for free ($) 



and periodic (H) boundary conditions one may show that 
there are two positive constants C($) and C(n) such that 



C($)|9A|<Tf) <0 



and 



Tr<T["^<C(n)|M|, 



(31) 



(32) 



where |i9A| is the area of the surface of A, i.e. the num- 
ber of terms in ^ which connect sites inside A to sites 
outside A. 

We have shown that a disordered quantum systems ful- 
fills a new correlation inequality which entails the same 
consequences as the first GKS inequality and holds in 
full generality without any restriction with respect to the 
classical case. It would be interesting to investigate cor- 
relation inequalities of type II (see H, lO] for the classi- 
cal case) as well as the validity of similar results on the 
Nishimori line 21| especially in view of the applications 
of the correlation inequalities to error correcting codes 
[2^ and their possible extension to the quantum case. 

The results we have presented in this letter can of 
course be extended to quantum Hamiltonian systems 
with general bounded interaction. 
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